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Abstract. We solve the recombination equation by tak- 
ing into account the induced recombinations and a phys- 
ical cut off in the hydrogen spectrum. The effective re- 
combination coefficient is parametrized as a function of 
temperature and free electron density and is about a fac- 
tor of four larger than without the induced recombination. 
This accelerates the last stage of the recombination pro- 
cesses and diminishes the residual ionization by a factor 
of about 2.6. The number and energy distribution of pho- 
tons issuing from the hydrogen recombination are calcu- 
lated. The distortion of the cosmic microwave background 
(CMB) spectrum depends strongly on the cosmological 
parameters f2, fiband-ffo and differs essentially from the 
Planck-spectrum for wavelengths < 190/im. 



Key words: cosmology: theory - cosmic microwave back- 
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1. Introduction 

The CMB radiation in a wide frequency range has a 
Planck-spectrum with T = 2.728 ± 0.002^ (Fixsen et 
al. 1996). The last dramatic event that could have in- 
fluenced a part of the spectrum was the recombination 
of the primeval hydrogen, since after it the residual ion- 
ization was very low and the radiation fields practically 
don't interact with the non-relativistic matter any more. 
The properties of the post-recombination Universe can 
be studied by observing the microwave background radia- 
tion. The fluctuations in the CMB are supposed to be the 
precursors of the largest structures observed today. 

To understand the nature of the microwave back- 
ground anisotropics it is necessary to have a correct pic- 
ture of the recombination process itself. The problem 
was first studied shortly after the discovery of the CMB 
by Peebles (1968) and at the same time by Zel'dovich, 
Kurt and Sunyaev (1968). Since than several authors have 
adressed the problem, but no one could improve on the 
basic approximations used in these works. However many 
details were worked out. Matsuda, Sato & Takeda (1971) 
examined the effects of the collision processes and found 



that they are negligible with respect to the radiation pro- 
cesses. Jones & Wyse (1985) improved on the calculations 
to allow for the presence of non-baryonic matter. Krolik 
(1989, 1990) showed that two previously equally neglected 
scattering effects in the Ly^ line almost completely cancel 
each other. Sasaki & Takahara (1993) concluded that tak- 
ing induced recombination into account changes the time 
history of the recombination at very low ionization grades, 
and lowers the residual ionization. They however used an 
approximation for the recombination cross sections and 
did not discuss the divergence of the sum of the recombi- 
nation rates. DeH'Antonio & Rybicki (1993) followed the 
progress of recombination numerically, using a multilevel 
hydrogen atom model; all angular momentum states were 
treated individually up to n = 10. Their recombination 
curve was similiar to that of Jones & Wyse (1985). In an 
other article Rybicki & Dell'Antonio (1994) calculated the 
time-dependence of the Ly a line profile in a homogenious 
expanding Universe. They found that the usual quasi- 
static approximation for the line shape is reasonable and 
does not cause substantial error in the solutions. They also 
determined the recombination history for several cosmolo- 
gies. Hu, Scott, Sugiyama & White (1995) discussed the 
effect of physical assumptions for CMB anisotropies and 
on the recombination process. 

Recombination, as it has been pointed out in the pa- 
pers quoted above, sets in when all energy levels of hy- 
drogen, except for the ground state are still in equilib- 
rium with radiation. When one hydrogen atom enters the 
ground state, a photon with energy greater than or equal 
to the energy of a Ly a transition is emitted. Photons, ca- 
pable of exciting electrons in the ground state, are lost 
via two main routes: 1. by cosmological redshift of Ly a 
photons, 2. by two-photon transition 2s — > Is. Of these 
two competing processes the first is purely cosmological, 
the second is atomic. Their relative importance depends 
on the cosmology. 

By theoretical examination of the recombination pro- 
cess one computes the time dependence of the ionization 
rate, the residual ionization, the position and width of the 
last scattering layer and the distortion of the CMB spec- 
trum. Though the amount of the residual ionization is 
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important for the further evolution of the Universe (Pee- 
bles 1993, Lepp & Shull 1984), it can not be measured 
directly. A possible observable consequence of the hydro- 
gen recombination in the early Universe is the distortion 
of the microwave background radiation spectrum. It was 
first calculated for the flat cosmological model by Peebles 
(1968). In this work we concentrate on the determina- 
tion of this distortion for different cosmological models. 
Recently there has been a considerable observational and 
theoretical activity to determine the spectrum of the cos- 
mic background radiation and carrying out galaxy counts 
in the far infrared/submillimeter range (Puget et al. 1996, 
Schlegel et al. 1997, Guideroni et al. 1997). Because the 
hydrogen recombination changes the CBM spectrum in 
the same spectral range a detailed exact recalculation of 
the frequency distribution of the recombination photons 
is important. In this paper we calculate the recombina- 
tion cross sections exactly by using a continuous, physical 
cut off for the highly excited states of hydrogen, take into 
account the induced recombination and explain, why the 
time dependence of the recombination process is so hardly 
effected by different techniques and by different effective 
ionization curves. 

The outline of the paper is as follows. In §2 we dis- 
cuss the recombination process and derive the recombina- 
tion equation. We give a new parametrization of the effec- 
tive recombination coefficients, taking into account the in- 
duced recombination as well. In the third section we solve 
the recombination equation. The spectrum of the recom- 
bination photons and its dependence on the cosmological 
parameters is given in the fourth §, and in the fifth closing 
section we discuss our results. 

2. The Recombination Process 

2.1. The Excited States of Hydrogen 

The number of recombinations in unit time can be cal- 
culated from the recombination equation as a function of 
the density of free electrons n e , temperature T , and cos- 
mological, and atomic constants. In this section we derive 
this equation following Peebles (1968). 

At the beginning of the hydrogen recombination the 
helium is already completely recombined. The mass frac- 
tion of the helium is 25 % of the total baryonic mass. As 
was remarked by Novikov and Zel'dovich (1967) the direct 
recombinations to the ground state are inhibited, while the 
new born energetic photons ionize again almost immedi- 
ately when there are already some hydrogen atoms. In the 
following we neglect completely the direct recombinations 
to the ground state. 

The states with principal quantum number n — 2 play 
a key role in the recombination processes. First we calcu- 
late the rate of recombinations to n > 2 for given free 
electron number density n e , temperature T and n 2s . Here 
n n i is the density of atoms in the state with the princi- 
pal quantum number n and angular momentum quantum 



number /. Second, we determine the net number of 2 — > 1 
transitions in unit time by given n 2s and n\ s . The two 
rates must be equal, so we can eliminate n 2s from the 
calculations. 

The binding energy of the n = 2 state is B 2 = Bi/4 « 
3.4 eV. During the recombination process there are a large 
number of photons with energy less than B 2l therefore the 
excited states of the hydrogen atoms are in thermodynam- 
ical equilibrium above the second level, i.e. 




«2 



= (2/ + l)e~^ rJS -n 2a . (1) 

In the case of gaseous nebulae the mean free path of low 
energetic photons is larger than the dimensions of the ion- 
ized region, the system is far from being in equilibrium 
above the n = 2 level in contrast to the recombining Uni- 
verse. 

The partition sum J2 nl (2l + 1 ) n (^ oltzmann ) j s diver- 
gent. As analysed by Hummer and Mihalas (1988) a num- 
ber of effects limit the range of summation. In the con- 
sidered temperature and density range the action of free 
protons turns out to be the most important, the ac- 
tion of neutral atoms are much smaller. The free protons 
destroy the state n with a probability 1 — w n , where 
w n = cxp[-(^ : ) 15 / 2 ] , with n* = 1075.9 cm" 2 / 5 n^ 2/15 
(n e in cm -3 ). This means that the recombined electrons 
become unbound with 1 — w n probability, before they 
would begin to move towards the state, corresponding to 
thermal equilibrium. The highly excited (n larger than 
~ n*) states are practically completely destroyed. The 
occupation numbers of the excited states in thermal equi- 
librium are: 

n nl = (21 + l)e" B2 "T n n 2s w n . (2) 

Owing to this fact, the number of electrons in bound states 
is finite, and thermal equilibrium between the n > 2 bound 
states and the continuum is possible. This approximation 
is good down to z w 300. The n = 2 levels freeze out 
between the redshifts 300 and 250. 

The ground state of the hydrogen atom is about 10.2 
eV deeper than the n = 2 level. Consequently, it's occu- 
pation is greater than the occupation of the excited states 
together: n\ s 3> X) n >2 n ™- This means, that the density 
of free electrons plus the density of hydrogen atoms in 
ground state can be taken as equal to the total proton 
number : n e + n\ s = p . The total baryon number (deter- 
mined by fib) is the sum of p and the number of baryons 
in the helium nuclei. 

2.2. Recombination on Excited States of Hydrogen 

The number of recombinations to the level nl and the 
number of ionizations from nl in unit time are 

dni n e n p = a n i nl and f3 n i n ni (3) 
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where a n i and f3 n i are the recombination and ionization 
coefficients, n p is the density of free protons. The usual 
definition of these quantities are 



I . (rcc) , Uon) 

a nl =< a nl V >Maxwcll, Prd =< O nl C >pi a nck 



(ion) 



(4) 



Here v and c are the velocities. The averages must be 
taken over the Maxwell-distribution for electrons and the 
Planck-distribution for photons. 
The ionization coefficients are: 



Pnl = {21 + 1) 



(ion) / 
C(T nl > 



8TT1S 2 



n3 „fe 



-dv . 



ew — 1 



(5) 



The relation between the recombination and ionization 
cross sections is: cr^j = <J^° n ^ -^r , where k and p are 
the photon and electron momenta. 

By using, instead of the electron velocity, the energy 
of the recombination photon in the Maxwell-distribution 
and expressing the recombination cross section through 
the ionization cross section, one arrives at 



. . . 4tt2(2Z+ f) s n 
c 2 {2um e kT) 2 



a^ ni (hvfe-^dhu,(6) 



B„ 



where m e is the electron mass. The integrals in a' nl and 
p n i are not the same. But, due to the principle of detailed 
balance (Milne f924) the recombinations and ionizations 
in thermodynamic equilibrium must exactly cancel each 
other. If we let contributions from the induced recombina- 
tions modify the recombination coefficients (e.g. Mihalas 
1984, Sasaki & Takahara 1993) instead of (§) we get 



Urd{T) = 



4tt2(2/ + 1) 

T e kT 

c 2 {2irm e kT) 2 

/ 2-Kh 2 \ 3/2 B- 



(ion) 
7 nl 



{hvf 



-d hi 



e w 



(7) 



The term — 1 under the integral in the denominator comes 
from the induced recombinations. Without it the principle 
of detailed balance would not be fulfilled. 

Summing up over all n, I levels and making use of equa- 
tion (J2J) the total rate of recombinations is: 



an,, 



with 



00 n—l 

EE 



Oinl W n . 



(8) 



The total rate of ionizations from the excited states of 
hydrogen is: 



00 n— 1 

EE 

n=2 ;=o 



Pni n n i w v 



kT\3/2 b 



. h 2 ) 



e zran 2s = /3n 2s . (9) 



The net number of recombinations to excited levels in unit 
time and unit volume is the difference of these two ex- 
pressions. The number density of the electrons changes, 
not only because of the recombination process, but also in 
consequence to the universal expansion. The simplest way 
to take this effect into account is to write on equation for 



the ionization rate, i.e. the ratio of free electrons to the 
number of free protons plus hydrogen atoms : 

d n fi 

- = a 

P 



dt p 



™k _ ( 



m P kT\ 3 / 2 b 2 T2, 2 
e kT 



2nfi z 



P 



(10) 



2.3. Depopulation of the n = 2 level 



Each 2p — > Is transition liberates a Ly Q photon, which 
can excite the ground state, and there are always enough 
low energy photons to ionize from an excited state. The 
transition to the ground state is final, when the number of 
these E y > B\ — B2 photons is also diminished. Photons 
can escape from the Ly a line either by redshift or by two- 
photon 2s — * Is transitions. 

In a stationary state the numbers of emitted and ab- 
sorbed Ly Q photons are equal. Consequently, the number 
of photons per mode in the Ly Q resonance line is: 



712 



n 2p 



3ni s — n 2p 3ni 



n 2p _ n^s 
nu 



(11) 



Because of the general expansion of the Universe a fre- 
quency 2/q is shifted in unit time by 



5vq 
St 



(12) 



where H is the Hubble-parameter. The number of recom- 
bination photons removed from the line by redshift in unit 
time is: 



Ris 



5v\ 2 %-nv\ 2 



St c 3 



(712 — 712 (background)) 



= H- 



'12 



c 3 O 12 e hu 12 /kT 



it- 



CIS) 



On the other hand, the rate of two-photon decay in 
the 2s states has been calculated by Spitzer & Greenstein 
(1951). The number of net decays in unit time is 



R 



2phot 



A, 



,(n 2s 



n ls e 



-{B 1 -B 2 )/kT\ 



(14) 



with A 2Si i s = 8.227s -1 . In thermal equilibrium the num- 
bers of two-photon decays and the two-photon Is — > 2s 
excitations are equal. 

At this point we have four unknown quantities: 
nis,n 2s , 712 andn e and four equations (|l0|), (|TT|) , ( |l3| ) and 
([Tij); so we are left with 



d n e 
dt p 



orn% 
P 



(3(1 



He 
P 



-(B 1 -B 2 )/kT 



%ttv1 2 c 3 H + A 2s ,i s (p- n e ) 
' 8irvf 2 c- 3 H + {A 2Sils + /3){p - n e ) ' 



(15) 



The time dependence of p, T and H are calculated from 
the Friedmann equation with the density parameter f2( = 
fib + f2 7 + ^ dm) j the baryon density parameter and 



the Hubble constant Hq - 
the cosmological constant) 



100hkms _1 Mpc -1 (A = 0, 
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2-4- Parametrization of the recombination coefficient 

To solve the equation ( |l5| ) we need the recombination 
coefficient a, wich is parametrized as a function of the 
temperature in several ways. Peebles (1968) used a = 
2.84-10" 11 /VTcm 3 s- 1 on the basis of Boardman's (1964) 
data for the seven lowest states. Zel'dovich, Kurt and Sun- 
yaev (1968) used 2.5 • Kr n /\/Tcm 3 s _1 . Peebles (1993) 
using Osterbrock's (1989) data got the parametrization: 
a = 4.1 • 10" 10 r-°- 8 cm 3 s- 1 . Recently Hummer (1994) 
tabulated very accurately the recombination rates. His re- 
sults confirm the parametrization given by Pequignot et 
al. (1991) in the form of 



a(T) = 1.261 • 10" 



T 



-0.6166 



1 + 5.087- 10- 3 T - 53 ' 



(16) 



Sasaki & Takahara (1993) had used the asymptotic form 
for the recombination cross sections which, also valid for 
large principal quantum numbers n, was less than 20% 
accurate for the important low lying levels. 

When calculating the recombination coefficients with- 
out taking into account the induced recombination, the 
— 1 falls out of the denominator of the integral in (0). 
In this case the a n series converges as log(n)/n 3 , so one 
can neglect the contribution of the terms with large n. If 
induced recombination is taken into account, the denom- 
inator under the integral in ((?]) goes to zero for large n, 
and the series of the a n diverges as 1/n . 

Here we calculate the total recombination coefficient 
from (||), wherein the factor w n , the probability that 
the states with principal quantum number n are not de- 
stroyed, makes the sum finite. The probability that a 
state is destroyed depends on the baryon density of the 
Universe. Consequently, the effective recombination coef- 
ficient also density dependent. 

Our procedure is then the following. We calculate the 
first several hundred a„z's using the exact results given by 
Storey & Hummer (1991). For large n (> 200), when the 
recombination coefficients have already reach their asymp- 
totic values (the difference between the exact and asymp- 
totic value < 0.1%), we use the a„'s given e.g. in Ry- 
bicki & Lightman (1979). The density dependence comes 
in through the w n probabilities. The effective recombina- 
tion coefficient is parametrized as a function of the tem- 
perature, T, and the electron number density, n e , in the 
form 



a(T,n e ) = 1.21 • 10~ 10 :zH 



3 o-l 



( 1 - 0.0226 hx(n e cm 6 ) ) cm 3 s 



(17) 



For a given temperature n e and a depend on the cosmo- 
logical parameters. However this dependence is mild, not 
larger than 10 % for the usual range of parameters. We 
compute the a(T,n e ) function together with n e by solv- 
ing (|l5|). The solution of this equation for n e and a(T, n e ) 
will be given in the next section. For the sake of clarity we 
present our recombination rate already here (Fig. 1). The 



form of the a(T, n e (T)) function is not far from the usual 
~ 7 l -0-7 function, but its amplitude is about four times 
larger than it would be without taking the induced recom- 
bination into account. Furthermore we compare our a(T) 
with (|l6|), and with Sasaki & Takahara (1993), whereby 
our function is somewhat steeper than that of these au- 
thors. 
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Fig. 1. Recombination coefficient as a function of redshift. 
Crosses: Present work (for Q — 1, fib — 0.1, h = 1); Continu- 
ous line: T~ ' 7 fit; Diamonds: Sasaki & Takahara 1993; Points: 
eq. (jl|); Dashed line: Peebles 1993. 



3. The solution of the recombination equation 

We solved the equation (|l^) for recombination numeri- 
cally using equations ([l7|) and (||). The results with and 
without taking into account the induced recombination 
are compared. 



3.1. The residual ionization 

At z ~ 200 H2 formation becomes possible. Since free 
electrons and protons serve as catalysts for the formation 
of molecular hydrogen (Peebles 1993), the value of the 
ionized fraction at this epoch is a very important. 

At about z ~ 500 (3 becomes much smaller than A2 B ,is- 
It is independent of the cosmology, because both [3 and 
A2 S is are atomic quantities. In this case the fraction 



in equation (15) is one. The second term in the square 
bracket is also small, so instead of the equation ([l5]) for 
x = ri e /p we can write 



ap 2 VU 2 

— = Jjq x = const— a; (1 



0.0226 hxpx). 



(18) 



dx 

dz zH flbh 

When the induced recombinations are not taken into ac- 
count the density dependent term does not appear and 
the residual ionization scales with VTitt^h^ 1 . The log- 
arithmic term contains another combination of the cos- 
mological parameters, therefore the scaling law is mildly 
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Table 1. Residual ionization at z = 220 for five cosmol- 
gical models. X2 and xi are the residual ionizations cal- 
culated with and without taking into account the induced 
recombination. 
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violated. The residual ionizations at z — 220 are shown 
in Table 1. When the induced recombination is taken into 
account the residual ionization is reduced by a factor of 
about 2.6 . Small deviations from this factor are due to 
the violation of the scaling law. The total number of ions 
is proportional to xSlb, so it is almost independent of the 
baryon density. 

We refined the previous treatment (Sasaki& Takahara 
1993) of the induced recombination by using the exact re- 
combination cross sections and applying a physical cutoff 
in the summation on the hydrogen states. This causes a 
10 - 20 % change in the residual ionization. 

3. 2. The time dependence of the fractional ionization 

The recombination history for different sets of comolog- 
ical parameters is shown in Fig. 2. The quotient of the 
fractional ionization at different cosmologies grows with 
time and approximately follows the x ~ \/fiQ^ 1 h~ 1 scal- 
ing law. We compare our result with DclP Antonio and 
Rybicki (1993), who did not take induced recombinations 
into account. For corresponding parameters our recombi- 
nation curve runs significantly under their curve, because 
the last phase of the recombination processes is acceler- 
ated by the induced recombinations. 

In Fig. 3 we show the ratio of fractional ionization cal- 
culated with (x2 ) and without (x\) taking induced recom- 
bination into account. In the first part of the process the 
induced free-bound transitions and the numerical value of 
the effective recombination coefficient have little effect on 
the course of the process, because the depopulation of the 
n = 2 states is much slower than the recombination to the 
excited states, the thermal equilibrium corresponding to 
the occupation of n = 2 is always maintained. Down to 
z w 1000 the x\/x2 ratio is one. Between z w 1000 and 



h = 1, £2= 1, £2 b = 1 

h = 1, £2= 1, £2 b = 0.1 

h = 1. £2= 1, £2 b = 0.01 

h = 1, £2 = 0.1, £2 b = 0.1 

h = 0.5, £2 = 0.3, £2 b = 0.01 
Dell' Antonio & Rybicki : 

h= 1. £2= 1, £2 b = 1 
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1600 1400 1200 1000 800 600 400 200 

cosmological redshift 

Fig. 2. Recombination history for different cosmological 
models. 
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Fig. 3. Ratio of the fractional ionization calculated with 
(£2) and without {x\) taking induced recombination into 
account. From right to left: 

a) = 1, 6 = 1, h = 1 

b) O = 0.1, 6 = 0.1, h=l 

c) = 1, 6 = 0.01, h = 1 

d) O = 0.3, 6 = 0.01, h = 0.5 



z ~ 500 the change owing to induced recombinations de- 
pends on the cosmology. Later, when z < 500, especially 
when the baryon density is small, the depopulation of 
n = 2 states is faster than the recombination to the excited 
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Table 2. Position and width of the last scattering surface, 
calculated for different sets of cosmological parameters. 
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Instead of the equation (|l^) we have 

Bx e x 



d d . . 



2ttHH(t)p(t) 



Q(x,t). 



(22) 



In the absence of sources the spectrum moves to the left 
(with "velocity" 1) in the (r , x) plane. 

When the right hand side of ( p2| ) differs from zero, the 
solution can be written in integral form as: 



states and this last process determines the net recombina- 
tion rate. At this stage the induced recombinations play 
an important role. As can be seen in Fig. 3, the x\Jx% 
ratio depends hardly on the cosmological parameters. 

3.3. The last scattering surface 

From the n e (z) function one can estimate the probability 
that a CMB photon has not been scattered since a given 
redshift z. Instead of this probability, usually its differen- 
tial is calculated. This quantity determines the probability 
density that the radiation was last scattered at z and is 
expressed as g(z) = e~ T dr/dz . Here r is the Thomson 
scattering optical depth: r = — J* co~Tn e {z)4^dz , and <7t 
is the Thomson cross section. The parameters, determin- 
ing the position and width of this layer are given in Table 
2. These results are in reasonable agreement with those of 
White, Scott & Silk (1994). This occurs because the last 
Thomson scattering happens early on, when the induced 
recombinations don't play any role. 

4. The spectrum of the recombination photons 

4-1. The time variation of the photon spectrum 

The evolution of the photon spectrum during recombi- 
nation can be calculated from the continuity equation in 
frequency space (Peebles 1968): 



d ( v n{v 1 1) 

dt \WT 



d 
ov 

v Qiy, t) 



vn(y, t) 



(19) 



n(u, t) is the number density of photons with frequency 
v, and Q{v, t) is the net rate of production of photons 
per unit volume and unit frequency interval. The Planck- 
spectrum fulfills this equation without the source term, so 
( |i~9] ) is also true for n'{v,t) = n(y,t) — Jipianck^, t) . Here 
we introduce new dependent variables with the definition 



5(M) = n'(v,t) 

and define the independent variables x, t as 



In- 



hv 
~B~i 



and 



H(t')dt'. 



(20) 



(21) 



9{x,t) 



x-\-t 



B x e x 



2irh,H(T>)p(T') 
with t' = x + t — x' 

1.2x10" 12 



-Q(x',T')dx', 



(23) 
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Fig. 6. Contribution of the Ly Q Photons and the 2s — > Is two-photon transitions to the distortion of the CMB for 
different cosmological models. 
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4-2. The spectrum of photons issuing from the 
recombination 

From the solution of equation ( |l5| ) we know the time 
dependence of the free electron density n e (t). Using the 
equations ([l4|) and (|l3| ) the net two-photon transition rate 
and the Ly a redshift rate can be calculated. The spectrum 
of the photons emerging from the 2s — * Is transition can 
be taken from Spitzer & Greenstein (1951). The appro- 
priately normalized two-photon spectrum in terms of the 
variable x is: 

$(x) = 0.7081e x ^/^e x \ (24) 

The ${x)dx = 2. 

For our purposes the frequency distribution of the Ly Q 
photons can be approximated by a delta function on the 
v axis, so 

Q(x, t) = R 2ph ${x) + R rs e x - Xl2 S(x - x l2 ) . (25) 

4-3. The dependence of the spectrum on cosmological 
parameters 

With n ls (t), the ratio of Ly Q redshift rate to the two- 
photon decay rate can be calculated from equations ( |l3|) 
and ([l4|). This ratio depends on the cosmology. For the 
2 — > 1 transition neglecting the — 1 in the denominator of 
equation ( |l3| ) and using the numerical values, one arrives 
at 

R _ R,s = &nvl 2 c- 3 H = 1 632.92 Q 1 / 2 
Riphot A 2s ,i s nu l-x (1 + zfl 2 Q b h ' 

The redshift of the new-born photons is given by equa- 
tion ( |22l ) in terms of the variables x and r. The spec- 
trum at different redshifts is shown in Figure 4 (fi = 0.1, 
Q b = 0.1, h = l) and Figure 5 (fl = 1, Q b = 0.01, h = 1). 
When the baryon density is low, the 2s — > Is transition 
plays a minor role. 

These results disagree with DeH'Antonio & Rybicki 
(1993), who state that the energy distribution of the pho- 
tons emitted by the 2s — > Is transition rather strongly 
peaked at v = rtLy Q /2 and the 2s — > Is transitions give 
no more than 1% difference in the free electron densities 
and the line strength. Our results (Fig. 4 and Fig. 5) show 
that the spectrum of 2s — > Is photons is broad (see also 
Spitzer & Greenstein 1951). Though for certain cosmolo- 
gies the contribution of these photons is small, but as it 
is demonstrated in Figure 6, for other set of parameters 
this contribution is quite considerable. In the models with 
high Q and low £l b the two-photon transitions are negli- 
gible, but for h 2 n b > 0.1 they are more important, and 
for the flat, high fl b model they are dominant. In spite of 
the above mentioned disagreement our recombination his- 
tory curve (Fig. 2) agrees well with that of Dell' Antonio & 
Rybicki (1993), because by computing this function they 
took into account also the 2s — > Is transitions. 



Figure 7 shows a part of the CMB radiation spectrum 
for different cosmological parameters, where the distor- 
tion due to recombinations is the most important. The 
hydrogen recombination begins when in the background 
radiation there are less energetic photons than hydrogen 
atoms. Consequently, the spectrum of the photons issuing 
from the recombination has a maximum near the Planck- 
curve. The other maximum in this spectrum corresponds 
to the two-photon transitions, is longwards from the first 
maximum at /ii/ max .2 = [B 2 — B±)/2 — 5.1 eV, and is 
under the Planck-curve. The 2s — > Is photons influences 
the distortion in two ways. First, the short wavelength 
part of their spectrum is above the Planck-curve and give 
a direct contribution to the distortion (Figure 6). Second, 
the number of recombinations are given by the number of 
hydrogen atoms. If a considerable part of recombinations 
happens by two-photon decays, there are less redshifted 
Ly Q photons, the amplitude and shape of their spectrum 
changes. 

For given Q b and h there are more photons above the 
Planck-curve when both fl and R are large. When fl and 
h are given, the spectrum with the larger lies above 
all the others. For a fixed value of and tt b the spectrum 
with larger h is larger. 

In Figure 7 we compare our spectrum with that of 
Dell' Antonio & Rybicki (1993). In the 120 /jm < A < 170 
/xm range there are about a factor of 2-3 more photons 
as seen by comparing our curve with the diamonds. That 
may be caused by the different technique in computing the 
recombination process. It seems us somewhat arbitrary to 
handel the 10 lowest level in a different way as the others. 
Burgess (1958) used this method to compute the recom- 
bination spectrum in nebulae. In the case of nebulae this 
method gives good results, because there is no thermal 
equilibrium above the n = 2 states, the electrons recom- 
bine on the low lying states and cascade down. 

5. Summary and Discussion 

Background photons stimulate recombination processes 
and that give an important contribution to the effective 
recombination coefficient. The resulting recombination co- 
efficient is about a factor of four larger than calculated 
without the induced recombinations (Figure 1) and de- 
pends, besides the temperature, on the free electron (pro- 
ton) density as well. In contrast to the recombining Uni- 
verse in a gaseous nebulae with the same temperature and 
free electron density the recombination coefficient is given 
by ([l6|), because there are no photons that should stimu- 
late recombination on excited states. The photons, except 
the Ly Q photons, emerging from recombination in the neb- 
ulae leave the nebula without interaction. 

For any given principal quantum number the a n i re- 
combination coefficients depend considerably on the angu- 
lar momentum quantum number I. There is a peak around 
I = 4 and a n i vanishes when / is high. Transitions between 
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states with high Ps are slow, which could influence the rel- 
ative occupation of s and p states. Since this last effect is 
probably small (see Hu et al. 1995), we have completely 
neglected it here. 

In the main part of the recombination process the free- 
bound transitions are more rapid than the depopulation 
of the n = 2 state and the thermodynamic equilibrium 
is maintained by reionization. Consequently, the course 
of the recombination process does not depend on the de- 
tails of the a(T,n e ) function. Because of that, authors, 
using quite different effective recombination coefficients, 
come up with similiar time dependence for recombination. 
When, however, the free electron density is already low, 



the depopulation of the n = 2 states is more rapid than the 
recombinations of the n > 2 states and the speed of the 
recombination process is determined by the free-bound 
transitions of the excited states and the induced recombi- 
nations are important. As a consequence, our ionization 
curve (Figure 2) is for small ionization grades steeper than 
that of DclP Antonio & Rybicki (1993), who did not take 
this effect into account. 

The last scattering of the CMB photons happens, with 
high probability, in the first phase of the recombination, 
when the number of free-bound transitions to the excited 
states is large. The induced recombinations and the nu- 
merical value of the effective recombination coefficient do 
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not influence the position of the last scattering layer. Our 
results agree with those of White, Scott & Silk (1994). 

On the other hand the residual ionization is a quan- 
tity determined by the details of the recombinations of 
the excited states at low n e . At z < 500, whether in- 
duced recombination is taken into account or not, the frac- 
tional ionization is proportional to VClCt^ 1 h -1 . However, 
because of the induced recombinations, the residual ion- 
ization is reduced by a factor of about 2.6. The exact value 
of this factor depends on the cosmological parameters. A 
consequence of the low residual ionization and the higher 
number of photons in the tail of the CMB spectrum is that 
the final abundance of molecular hydrogen will be much 
smaller than that estimated by Lepp & Shull (1983). 

The hydrogen recombination gives contribution to the 
CMB spectrum in the A < 190 /im range. The amplitude 
and form of this distortion depends on the number of pho- 
tons issuing from the recombination, and their distribution 
between the two decay way, the 2s — > Is two-photon mode 
and the redshiftcd Ly a mode. The induced recombinations 
only in the last phase influence the recombination process, 
when the ionization grade already very low. This gives a 
very small contribution to the distortion of the spectrum. 

The contribution of helium recombination to the pho- 
ton spectrum was discussed by Lyubarsky & Sunyaev 
(1983) and by Fahr & Loch (1991). According to their re- 
sults the recombination of He sets in at a redshift of about 
four times larger than that for hydrogen recombination. 
But the photon energy is also four times larger and the 
photons emerging from helium recombination could also 
distort the CMB spectrum in the same frequency range as 
the H-recombination. However, as it was pointed out by 
Peebles (1995), at those redshifts where He II recombines 
there is already a trace of recombincd hydrogen, which 
can absorb some of the photons created during helium re- 
combination. Thus direct recombination into the ground 
state is possible for helium and the Saha formula is a good 
approximation to describe the helium recombination his- 
tory. The optical depth due to ionization between zi and 
z 2 is 

f Z2 dr 

T = J (Tionnu—dz. (27) 

The redshift of the helium recombination is zi, 
the cross section for hydrogen ionization is <Ji on [rts 
(Bi/hv) 8 ^ 3 , with v(z) = v(zi)j^}, and the number density 
of hydrogen atoms, as computed from the Saha formula, 
is ni s (z) . The distance to the radiation source at the mo- 
ment of emission, r(zi), is calculated from the Fricdmann 
equation. The integral reaches 1 at about z 2 ~ 1900. Con- 
sequently, in contrast to the results mentioned above, he- 
lium recombination does not leave any trace on the CMB 
spectrum. 

As is discussed in Hu et al. (1995), the difference be- 
tween the temperature of the photons and the kinetic tem- 
perature of the electrons has an effect below z 200. We 



follow the recombination process only down to z = 220, 
so this effect is neglected. 

The possible ns — * Is and the nd — > Is two-photon 
transitions could also have some importance. The number 
of bound electrons in the n > 3 states is of the same order 
of magnitude that for n2 S . The two-photon decay proba- 
bilities for these states are smaller than for ri2 S • Moreover 
some of these photons will have more energy than B2 - Bi 
= 10.2 eV. However, the contribution of these transitions 
is not entirely negligible, they can change the balance be- 
tween the two routes, redshift and two-photon decay, for 
the elimination of energetic photons and change the spec- 
trum of the recombination photons. This correction will 
be discussed in detail elsewhere. 

If there were some energy input prior to the hydro- 
gen recombination and a " y-distortion" of the spectrum 
(Zcl'dovich & Sunyaev 1969) it would be easier to observe 
the consequences of H-recombination (Lyubarsky & Sun- 
yaev 1983). However, the observations permit a very low 
upper limit for the y-distortion (Mather et al. 1994, Fixsen 
et al. 1996). 

The distortion seems to occur at wavelengths where 
other sources give considerable contributions. By doing 
the calculations we had hoped that the zodiacal forground 
emission, the dust and molecular emission from the in- 
terstellar medium would be reliably modeled, subtracted 
and the calculated spectral distortion would be measur- 
able. Recently that substraction has been done (Schlegel, 
Finkbeiner & Davis 1997) and it turned out that in the 
considered wavelength range the measured upper limit of 
the background is at least two order of magnitude higher 
as the spectrum calculated here. As it was pointed out 
by Guideroni et al. (1997), it could originate from early 
starlight scattered by dust. 
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